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We calculate the angular momentum and energy of a vortex dipole in a trapped atomic Bose- 
Einstein condensate. Fully analytic expressions are obtained. We apply the results to understand a 
novel phenomenon in the MIT group experiment Q], an excellent agreement is achieved, and further 
experimental investigation is proposed to confirm this vortex dipole mechanism. We then suggest 
an effective generation and detection of vortex dipole for experimental realization. Application of 
the sum rule to calculate collective mode frequency splitting is also discussed. 

PACS number: 



I. INTRODUCTION 

The distinguished behavior of superfluid to carry an- 
gular momentum have attracted people's attention from 
long ago. It is well known that vortex excitation plays 
an important role when a superfluid responds to the ex- 
ternal rotation. Nowadays, the vigorous development in 
cold atom physics has provided a new and more favorable 
arena to prove many theoretical predictions of rotating 
superfluid's behavior. It has been clearly demonstrated 
by many recent experiments that the vortex will be a 
stable excitation with quite long life-time when sufficient 
angular momentum is imparted, and the vortex lattice 
will be generated when even more angular momentum is 
brought to this system. 

Besides the vortex itself, the vortex dipole, i.e., a 
bound vortex-antivortex pair, is also of particular in- 
terest. It is well known that the low- lying behavior of 
two-dimensional helium superfluid is dominated by vor- 
tex dipole, for the excitation energy of a vortex dipole 
is much smaller than a single vortex. With the increase 
of temperature, the vortex pairs will become unbounded 
and vortices will be free. This will lead to a phase tran- 
sition known as Kosterlize-Thouless(KT) transition. Al- 
though no vortex-antivortex pair has yet been observed 
directly so far in atomic Bose-Einstein condensate(BEC), 
it is natural to presume that vortex dipole may play an 
important role when the imparted angular momentum is 
not large enough to support the vortex excitation. So 
it is worthwhile to study the properties of vortex dipole 
excitation in an atomic BEC. 

The existence and importance of vortex dipole in an 
atomic BEC are also indicated by some recent experi- 
ments and numerical simulations. MIT group stirred the 
BEC locally with the stirring frequency which was too 
small to generate any surface mode, however, they ob- 
served a large number of vortices after stirring for some 
time. This suggested a local mechanism of generating 
vortex, say through vortex dipole excitation[l|. Recently 
some numerical calculation also indicated that stable vor- 
tex dipole could exist in a trapped atom BEC for a broad 
range of parameters |3- All of these lead us to believe 
that cold atom physics will offer people opportunity to 
investigate this intricate vortex dipole excitation in su- 



perfluid. What is more, the characteristics of trapped 
atomic BEC will add much richer physics content to the 
properties of vortex dipoles, as well as to the experimen- 
tal detection and production of vortex dipoles, which will 
be emphasized in this paper. 

To proceed in this interesting and meaningful research 
field, we of course should first make clear the properties 
of a single vortex dipole in BEC, i.e., its angular mo- 
mentum and energy. In section II, we will first deduce 
an exact result of its angular momentum and energy. It 
is not trivial to calculate the angular momentum and 
energy of vortex dipole in an inhomogeneous condensate. 
Here we use an integration method in complex coordinate 
to obtain fully analytical results, which will be described 
explicitly in this section and the appendix. We then ap- 
ply our results to investigate MIT group experiment 
Through our calculations, we provide strong evidence to 
support the presence of vortex dipole during this kind 
of vortex nucleation process, and suggest further experi- 
ment to confirm this mechanism. 

For experimental study the properties of vortex dipole 
in details, an effective way to generate and detect it is 
required. Although various methods have been widely 
adopted to produce vortices successfully, the formation 
of vortex dipole has not been systematically discussed. 
In section III, we will suggest that JILA's method 
which generate vortex through a coherent spatial and 
temporal coupling between two-component BEC, could 
also be used to produce vortex dipole. 

Measuring the density distribution could demonstrate 
the density depression in the neighborhood of the vortex 
center but have nothing to do with its phase. The same 
problem will be encountered in the case of vortex dipole. 
We could not distinguish which vortex is positive or neg- 
ative only from the density image. What is more, the 
velocity field of a vortex dipole decreases quickly from 
the dipole center, and any closed path, which is a little 
far from the vortex dipole and includes both the vortex 
and antivortex, will be quite similar with the case that 
there is no singularity at all in this area. It is also not so 
effective to use interference method to reveal the phase 
information of vortex dipole as in vortex caseQ. How- 
ever, experimentalists now could accurately measure the 
frequency splitting of the collective modes when the time 
reversal symmetry is broken due to the presence of vor- 
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tices. This splitting will clearly distinguish vortex dipole 
from single vortex and vortex-free states. In theory, there 
are usually two ways to calculate the frequency of the col- 
lective mode in presence of vortices, one is perturbative 
calculation and the other is the sum rule approach. The 
results obtained from these two approaches are in good 
agreement with each other when the vortex is in the cen- 
ter of condensate. However, they are no more consistent 
when the vortex is off-center as pointed out in Ref.|j|. 
In section IV, we first use perturbation approximation to 
calculate the frequency splitting of the collective mode, 
and then we explain why the sum rule approach fails here. 
A brief conclusion and discussion will be given in the last 
section. 



II. PROPERTIES OF A VORTEX DIPOLE 



Only in the region r\ < r < r 2 , 9 will change 2n after 
going along a circle. Thus, 



L = 2ith 



prdr. 



(6) 



we notice that the energy of a vortex dipole increases with 
their distance, and the most efficient way to carry angular 
momentum is that f\ and f 2 direct along the same ra- 
dius, namely they are fully polarized, as are schematically 
shown in figQ). So we could denote r\ = D — d/2,r 2 = 
D + d/2, therefore 



Dd D 3 d d 3 D 



L = 2Kh Po R 2 {—-—-- 
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(7) 



where po is ftio^mR 2 /(2g). Now we turn to the energy of 



A. Angular Momentum and Energy 

We hereafter consider a BEC confined in a pancake 
harmonic trap, which allows us to deal with a quasi-two- 
dimensional system. The Gross-Pitaevski equation in a 
non-rotating coordinate frame is 

_^vV + ^ + slV'l 2 ^ = ^, (i) 

Zm 



where V — \mio\r 1 , and g is the effective quasi- two- 
dimensional interaction strength, obtained by integrating 
the usual three dimensional one Airh 2 a sc /m over the z di- 
rection. As we are considering the small angular momen- 
tum case, we could start with the Thomas-Fermi(T-F) 
approximation. So the density distribution of the ground 
state is well known as 



, . Hlu 2 , m , n 9n 
P{r) = ^—(R 2 ~r 2 , 



(2) 



where R is the T-F radius. Once a vortex dipole appears, 
in the region outside of vortex core we can write 



with 



9(r, tp) = Q(f— r*i) — 0(r — r%). 



(3) 



(4) 



Here O(r) is the angle of r with respect to x axis. r% 
and r 2 denote the positions of vortex and antivortex re- 
spectively. Eq.lPJ follows from the assumption that the 
presence of vortices will not change the density distribu- 
tion of the condensate but only introduce a phase, except 
inside the vortex core. This is true when the vortex core 
A is small in the T-F region. With Eq.© we could cal- 
culate the angular momentum of the condensate as 



L = h 



prdr 



dp 



dp = H prdr6(r,tp)\f.(5) 




FIG. 1: Vortex dipole configuration: the grey torus expresses 
the effective part of the condensate which carries the angular 
momentum 



the vortex dipole. In the T-F region, it could be written 
as 

h 2 



where 



E = — / pv dxdy, 



z x (r — z x (r — ?2) 



\T — r\\ 



\r - r 2 \ 



(8) 



(9) 



Here we have neglected the presence of mirror dipole due 
to the edge effect, which would be a good approximation 
in a broad region except when the dipole is very close to 
the boundary of the condensate. For convenience, we do 
all the calculation in complex coordinates. We will use 
z — x + iy and z = x — iy, and then obtain 



h 2 p 



(R 2 -\z\ 2 ) 



\zi - z 2 \ 



AimR? jq 
According to the Green formula 



\z - zi\ 2 \z - z 2 \ 2 



<)!> 



dzAdz. (10) 



(11) 



we could obtain the energy expression in the region A <C 
d< R 



E = 



2h 2 p 



In- 



A' 



(12) 
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The detailed calculations are given in Appendix B. As the 
energy is dependent on d logarithmically, the vortex and 
antivortex should bind up to each other closely, and we 
could neglect higher orders of d/R. Compared with the 
well known result for homogenous case, we also find that 
in Ea. l)12|l the homogenous density is replaced by local 
density at the position of vortex dipole. It can be easily 
understood because the amplitude of the velocity field of 
vortex dipole falls rapidly from the center of the vortex 
dipole, and most contribution to the energy of a dipole 
comes from the local velocity field. 



B. Interpretation of Experiment 

It is of long standing interest to investigate the intrinsic 
mechanism of vortex nucleation since the study of helium 
superfluid. In Ref.Q dynamic nucleation of vortices in a 
trapped atomic BEC was studied experimentally. When 
the condensate was stirred by a laser with the beam waist 
comparable to the T-F radius, the condensate was glob- 
ally rotated. In this case they found an enhanced vortex 
generation when the stirring frequency coincides that of 
surface excitations. It is indicated that vortex could be 
generated through surface excitation. However, when the 
condensate was stirred locally by a small beam, and the 
stirring frequency was far below any surface excitation 
frequency, they also found large number of vortices. This 
phenomenon suggested that the vortex can also be cre- 
ated locally in the bulk of the condensate, which seems to 
be in conflict with the topological argument, unless the 
vortex dipole was excited as an intermediate step. With 
the increase of angular momentum, the anti- vortices will 
finally be expelled out from the system and leave only 
the vortices in the condensate. 

A remarkable feature of the local nucleation is that the 
vortex number produced by this local stirring is strongly 
dependent on the distance I between stirring position and 
the center of condensate. The vortex number presents 
a maximum in the region where the ratio of I to the 
condensate size R is around 0.4. Hence, a question can 
be raised that whether we can understand the maximum 
from the properties of vortex dipole excitation. 

Our answer to this question is, it is just in this region 
that the excitation energy of vortex dipole takes a min- 
imum. In other word, the vortex dipole is easiest to be 
excited in this region. With Eq.Q and Ea. l|12|) in hand, 
we could calculate the energy as a function of the posi- 
tion of vortex dipole D for a given small value of angular 
momentum Lq- The result is 
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2nhpoB? A(D) R - (| )3 



(13) 

where Eq denoting 2h 2 po/m is introduced as an energy 
unit. Here we should emphasize that the vortex core size 



A also depends on D/R via 
A(D) = 



Ar 



l-(f) 2 



(14) 



This result is plotted in Fig©. We could clearly see 
that the energy minimum occurs at D/R around 0.4. Our 
calculation thus provides a support to the vortex dipole 
mechanism. 
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FIG. 2: Vortex dipole energy. Solid, dashed, and dash dot de- 
note T = 0.05, 0.1, 0.2 respectively, where T = L /2-Khp RK 



To further confirm the vortex dipole mechanism, we 
investigate the properties of vortex dipole in an anhar- 
monic trap 



V(r) 



1 2 2,1 , 2 4 

—rauj i r H — mk r . 

2 ^ 2 



(15) 



The angular momentum and energy of the vortex dipole 
will become 



D 5 d 
~R?R 




(16) 



m x . (ir) 



where A = i . We find that with the increase of the 

quartic term, the minimum in the E — D curve moves 
toward larger value of D/R as shown in fig.® Hence we 
predict that the stirring position where the number of 
generated vortices reaches its maximum will also change 
correspondingly, if a similar experiment is performed in 
an anharmonic trap. 

III. VORTEX DIPOLE FORMATION 

Let us first briefly review the vortex formation 
by phase imprinting method^, which used a laser 
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FIG. 3: Vortex dipole energy. Solid, dashed, and dash dot 
denote A = 0.5, 1, 2 respectively, and T is fixed at 0.05 

to couple a two-component BEC. When we write 
the order parameter ijj as n 1 / 2 e m %, where x T — 
(cos(6»/2)e-^/ 2 , sin(6»/2)e^/ 2 ), we could map the wave 
function to a pseudo-spin % on a Bloch sphere. Under 
the optical coupling, the pseudo-spin will rotate with the 
effective Rabi oscillation frequency fl e s- In the experi- 
ment the laser was locally applied and moved around a 
circle with the frequency uj r , and the condensate in dif- 
ferent points acquires a phase difference proportional to 
ficffAi, where At is the time interval of the laser pass- 
ing through these two points. If uj r = fl e g, the phase 
difference between the two points is identical to AO, and 
the phase configuration of a vortex is built up. Besides, 
if we only change the sign of O e g and keep the circling 
direction of the laser unchanged, an anti-vortex will be 
produced. 

Extending the above method to dipole generation is 
straightforward. From Fig.QJ, we can see that the ve- 
locity field of a vortex dipole is clockwise in left plane 
and anti-clockwise in right plane. To produce a vortex 
dipole, the moving laser is required to change the sign of 
f2 e ff once it crosses the y axis. What is more, the phase 
configuration of the vortex dipole could be produced ex- 
actly. f2 e ff, the velocity of the moving laser v r and the 
phase change gradient of the vortex dipole V# should 
satisfy the relationship 

v r ■ V0 = n cS (is) 



everywhere along the circle. 8 denotes dipole phase ex- 
pression Eq. Q . Once the phase along the circle is fixed, 
according to the uniqueness of solution to the equation 
V 2 (9 = 0, the two phase singularities, as well as the phase 
elsewhere inside the circle will consequently be deter- 
mined. 
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FIG. 4: Vortex dipole generation. Left: arrowhead denotes 
the velocity field of the vortex dipole along the circle. Right :v r 
denotes the velocity of the moving laser 



IV. VORTEX DIPOLE DETECTION 

When the ground state of the condensate possesses 
time reversal symmetry, the ±m collective modes will be 
degenerate. However, once the time reversal symmetry 
is broken due to the presence of vortex, the frequencies 
of the two modes will split. Zambclli and Stringari used 
the sum rule approach to obtain an analytical result of 
this splitting when a single vortex lies in the center of the 
condensate_6| . Svidzinsky and Fetter also obtained a con- 
sistent result by solving the the hydrodynamic equations 
with perturbational approximation^- They both found 
that the frequency splitting is proportional to the total 
angular momentum carried by the condensate. Later the 
frequency splitting when the vortex is off-center was also 
calculated pertubatively 8] 5] . They found that the result 
obtained from the sum rule approach does not agree with 
that from the perturbative approach any more. In this 
section, we will first calculate the collective mode split- 
ting in presence of a vortex dipole from the perturbation 
approximation, and we find the similar disagreement be- 
tween perturbation approach and the sum rule approach 
will also occur here. We will briefly discuss the reason 
why the sum rule approach fails to meet the result of per- 
turbative approach. We should point out that it is not 
the sum rule itself but the incorrect way we use it that 
leads to this failure. 



A. Frequency splitting induced by vortex dipole 

Substituting ip — ^fpe 1 ^ into Eq.JD, we could obtain 
the well known hydrodynamic equations in the T-F re- 
gion. 

^ + v-(H = o, (is) 

Q — * -j 

M-^+V(V + gp^ f i+-Mv 2 ) =0, (20) 
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where v = V0. Expanding p = po + 5p and v = vq + 5iJ 
with 

zx (f-fi) z x (r — ro) 
v = vi + v 2 = ^,9 - ,^ , (21) 



r - n 



r - r 2 



and after linearization, we could obtain the equation for 
collective modes 



(Lu 2 +-^V-p V)6p 



-iuf{v 5p)+^\/-[p V(<ij -f)6p]. 

ML) 



(22) 

We notice that the left hand of Ea. 1)22(1 is just the equa- 
tion of collective motion corresponding to a vortex-free 
ground state. And we also note that the effect of vo to 
the collective modes is always small, because vo decreases 
quickly from the vortex dipole center as 1/r 2 while Sp 
increases from the condensation center as r m , their over- 
lapping is always small except the position of the vortex 
dipole is nearby the boundary. So we could treat the 
right hand side as perturbation safely and obtain 



lo 2 - = -iu> (g\v ■ V\g) + 



1-9 



(g\V-p V{v -V)\g), 
(23) 

where loq is the collective mode frequency in a vortex free 
case. Hereafter we consider the following surface mode 



Iff) 



1 r" 



R m+1 



(24) 



From the calculation presented in Appendix C, we could 
obtain the splitting between +m mode and —m mode as 



LU + — D- 



2(m + l) 
IT 2 



2m -2 



(25) 



In experiments, from the optical absorption image the 
vortices positions r\ and r% can be determined, and then 
by measuring the frequency splitting and comparing the 
obtained result with Ea. (|25l) . one can distinguish a vor- 
tex dipole from two vortices with the same winding di- 
rections. 



B. Discussion on the Application of Sum Rule 

When we apply the formulism of Ref.[6j to calculate 
the frequency splitting of ±m modes, we can not obtain 
the same result as Ea. (|25|l . To understand the reasons, 
let us first briefly remind ourself of the sum rules ap- 
proach. The general formula of the sum rules for an 
arbitrary operator T can be easily proved as 



(26) 



If we choose k = 0, and with a particular choice of the 
operator T which satisfies |.F n o| ~ S u q, we could obtain 



fc2 2 

h w„ = 



[[F+,H],[H, M|) 
{0\[F+,[H,?\]\0) 



(28) 



This is a powerful method because we can obtain the ex- 
citation energy directly from the property of the ground 
state without the explicit wavefunction of excitation 
state. However, the successful application of the sum 
rule approach relies on the proper choice of T, which 
should excite the ground state to only one definitive ex- 
cited state. When the system possesses some particular 
symmetry, we may easily find out a suitable T . That's 
true when there is a central vortex, as rotation invariance 
is still present and the excited state should also have the 
form e lm ^. So we could easily select T as (x ± iy) 2 to 
make use of the sum rule. However, the rotation invari- 
ance is broken in the case of vortex dipole. The choice of 
T still as [x ± iy) 2 will excite several states, and Ea. p%|) 
will not tenable. That's why the sum rule approach fails 
to agree with the perturbative result in Ref.[5j. 



V. CONCLUSION 

In summary, by applying a convenient integration 
method in complex coordinate, we have obtained the 
fully analytic expressions of the angular momentum and 
energy of a vortex dipole in a trapped BEC. We also 
suggest an effective method to generate and detect the 
vortex dipole. Our work may provide a starting point for 
the research in the problem of vortex dipole excitation. 
There still remains many interesting fundamental issues 
in both theoretical and experimental physics, such as its 
dynamics, stability and detailed generation mechanism. 

The precession of off-center vortex in a trapped con- 
densate has been fully investigated both in experiment 
and theory. The mutual interaction between vortex and 
antivortex is presumed to bring a more nontrivial tra- 
jectory for vortex dipole. On the other hand, its dy- 
namic behavior is closely related to its stability. In the 
dual electromagnetic picture t ne positive and nega- 
tive charge, to which the vortex dipole corresponds, will 
tend to annihilate each other via electromagnetic radia- 
tion, say phonon excitation in original picture. However, 
in a trapped condensate, there are two mechanism to 
stabilize the vortex dipole, one is the angular momen- 
tum conservation, and the other is the discrete spectrum 
of the the phonon in a finite system. Investigation of this 
interesting competition will reveal the underlying physics 
of the numerical simulation^. The fundamental mech- 
anism of generating vortex dipole from a local stirring 
is even more complex, and it may resemble the birth of 
electron and positron from the vacuum polarization of 
the electromagnetic field. 

Another important issue which should also arouse 
great attention is the complex phase diagram of rotating 
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two-dimensional superfluid. When the angular momen- 
tum L equals zero, the critical temperature T c for K-T 
transition is well understood. But when the superfluid 
carries a certain angular momentum, the dependence of 
T c to the angular momentum, as far as we are concerned, 
hasn't been clearly revealed. On the other hand, with 
the increase of angular momentum, the vortex dipoles 
will break into free vortices, and finally stable vortex ex- 
citation will become dominative. However, the explicit 
mechanism of the transition from vortex dipole excita- 
tion to stable vortex excitation is also still unknown. 
The interplay of the angular momentum and the ther- 
mal excitation will bring rich physics to the rotating two- 
dimensional superfluid. Several works are proceeding. 
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APPENDIX A: INTEGRATION IN COMPLEX 
COORDINATE AND GREEN FORMULA 



APPENDIX B: VORTEX DIPOLE ENERGY 

To obtain the result of the integration Eq. (|10f> we first 
calculate 



\zi - z 2 \ 



Q \z — Z\\ \z — Z 2 \ 2 



dz A dz. 



According to above method, it turns to 



zi - z 2 



In -dz 



(z-z 1 )(z-z 2 ) z-z 2 



(Bl) 



(B2) 



We notice that not only does the integrand has two 
singularities z\ and z 2 , but also includes multi- value 
parts ln(z — z\) and ln(z — z 2 ), so the integration in the 
boundary should be as the left side of Fig-.J^J. 




FIG. 5: The integration along the boundary in the complex 
coordinates. 



In this section we introduce our method used to do 
integration in a two-dimensional space. By denoting z = 
x + iy and z = x — iy, we can turn the integration into 
the complex coordinate as 



/ f(x,y)dxdy = f /( 
Jn 2i Jn 



z, z)dz A dz, 



(Al) 



where 51 denotes the integration region, and the A op- 
erate fulfils dxi A dxj = (Sij — l)dxj A dxi, here Xi is 
the coordinate in real or complex space. Introducing a 
differential 1-form uj = g(z, z)dz in which the function g 
satisfies 



dgjz, z) 
dz 



f(z,z), 



(A2) 



we have duj — f(z, z)dz A dz. According to the Green 
formula, we could change the above integration over 
to over its boundary denoted as Oil as 



f(z, z)dz A dz 



g(z,z)dz 



(A3) 



i){>. 



This method is very useful when doing integration in a 
irregular region, because it turns a two-dimensional in- 
tegration into a one-dimensional one, and result of the 
latter is easier to be obtained numerically or approxima- 
tively. 



In the condition that d ^ A, the integration along 
circle C\ reads 



z\ - z 2 



In -dz 



5l =Ae«o ( z ~ z i)( z ~ z i) z ~ z i 

d -— -[In ^ - i(ip + n)]i(z - z^dip 



= — 2«7rln 



A 



The similar result of integration along C 2 circle reads 
2in In 4 . We notice that the imaginary part of the log- 
arithmic function in integrand in branch l\ differs from 
that in l 2 by 2-7T, the integration along the two branches 
gives 



d-A 



-d „ , , d 

—zirdx — 4«7rln — . 

A x(x — d) A 



(B3) 



The integration is more complex along circle C3, we shall 
decompose Ea. (|B2|) as 



i=l,2 



ln(z - n) 



z — n 



dz 



\a{z — rj ) 



dz. 



(B4) 



The general form of the integrand could be written as 



i + E(|) r 

m—l 



In z 



E(fr 



(B5) 
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then each of the four terms after decomposition is easier 
to be calculated. We obtain the integration along C3 
circle as 



2in 



In 



1 



R 2 



1 



-2 
R 2 



1 - 



\ziz 2 \ 
R 2 



(B6) 



The above value will vanish in the limit of djR — > 0, thus 
could be neglected. Then we shall obtain 



\zi - z 2 \ 



T;dz A dz — 8i7r In — . 
2 A 



In \z — z\\ \z— zi\ 
Proceeding the similar steps, we could obtain 

r \ z \ 2 \zi-z 2 \ 2 

Jn 



, , D\ d 

7?dz Adz = In — . 

2 R 2 A 



in R 2 \z — z\\ \z — z 2 
The energy of a vortex dipole will finally read 



E = 



2h 2 Po 
M 



1 d 



(B7) 



(B8) 



(B9) 



APPENDIX C: FREQUENCY SPLITTING 

We notice that the Eq. l|23[l is linear with respect to 
vo, we could first calculate a single positive vortex case. 
That is 



2iz'z' ' 



and 



(CI) 



(C2) 



11 2z>z> ' 

where z' — z — z\, z\ denotes the location of vortex. Thus 
(9\vo-V\g) = {g\ T ^ ^-\g) 

z — ZZ\ z — Z\ 



m + 1 
\R 2m + 2 

m + 1 
\R 2m + 2 



n z - z i 



-d z z m dz A dz 



-dz 



an z ~ z i 



(C3) 



The integration along the boundary is expressed as the 
right ride of FigJSJ), and the result is easily to be obtained 
as 



G?|wo • V\g) = 



(m + l)i 
~~R 2 



1 - 



R 2r 



(C4) 



With the same process, we could also obtain (g|V 
/9 V(u • V)\g) as 



2im(m + \)po 
R 1 



1 



l£i| 
R 



2m 



- 2 



|£i| 
R 



2(m-l)' 



(C5) 

Recalling that u>q — u>±y/rn and substituting the above 
expression back into Ea. (|23|l . we will obtain 



LU — LUq = 



m + 1 



1 - 



|£i| 
R 



2(m-l)' 



(C6) 



The same result of Ref.( 5]) is again obtained. As we 
have emphasized that Ea. l|23|l is linear with vq, the above 
result of single vortex can be directly extended to vortex 
dipole, that is 



bJ — UJq 



m + 1 



!£2j 
R 



2(m-l) 



ill 

R 



2(m-l) 



where z 2 is the position of anti-vortex. To calculate the 
frequency of — m mode, we only need to change the sign 
of vq equivalently. So the frequency splitting of ±771 mode 
is 



UJ + — CJ_ = 



2(m+l) 
~R 2 



l£2| 

R 



2(m-l) 



l£i| 
R 



2(m-l) 



(C8) 
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